set more off
estsimp logit forcea civila civilb  insta instb relpow jointcont jointally jointdem  instb_civilb insta_civila peaceyrs

setx civila 1 civilb 0  insta 1 instb 0 instb_civilb 0 insta_civila 1  relpow mean jointcont 1 jointally 0 jointdem 0 peaceyrs mean, 
simqi, prval(1) genpr(civilwarinstabhigh)
setx civila 0 civilb 0  insta 1 instb 0 instb_civilb 0 insta_civila 0 relpow mean jointcont 1 jointally 0 jointdem 0 peaceyrs mean
simqi, prval(1) genpr(instabhigh)
setx civila 1 civilb 0  insta 0 instb 0 instb_civilb 0 insta_civila 0  relpow mean jointcont 1 jointally 0 jointdem 0 peaceyrs mean
simqi, prval(1) genpr(civilwarhigh)
setx civila 0 civilb 0  insta 0 instb 0 instb_civilb 0 insta_civila 0 relpow mean jointcont 1 jointally 0 jointdem 0 peaceyrs mean, 
simqi, prval(1) genpr(civilwarinstablow)

gen fd_1=civilwarinstabhigh-instabhigh
gen fd_2=civilwarhigh-civilwarinstablow

centile fd_1, centile(2.5 97.5)

centile fd_2, centile(2.5 97.5)

gen sec_fd=fd_1-fd_2
centile sec_fd, centile(2.5 97.5)

drop b*
